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Abstract 



In this paper, we report on the optical absorption in porous 
sihcon. We model the absorption process assuming that porous silicon is a 
pseudo ID material system having a distribution of band gaps. We show 
that in order to explain the absorption we specifically need to invoke - (a) 
k is not conserved in optical transitions, (b) the oscillator strength of these 
transitions depend on the size of the nanostructure in which absorption takes 
place and (c) the distribution of band gaps significantly influences the optical 
absorption. We also show that it is not possible to extract the band gap of 
porous silicon from a plot of y/ahuj vs hu). 
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I. INTRODUCTION 

The discovery of efficient photoluminescence^^'"^^ from porous silicon (PS) has 
attracted the attention of many researchers. The photoluminesecence peak occurs at an 
energy greater than the band gap of c-sihcon and can be tuned through the visible 
spectrum by changing the preparation conditions. Canham'^^^ proposed that PS is a 
nanostructured material and the band gap of these nanostructures is enhanced due to 
quantum size effects. He argued that this would account for the fact that the luminescence 
energy is greater than the band gap of c-silicon. Evidence for the increase in band gap was 
first reported by Lehmann - from optical absorption measurements*-^-* . Transmission 
electron microscopy measurements reveal that the diameter of the columnar 
nanostructures that make up PS depends on preparation conditions and can be as small as 
2-4 nm*^^'^^. PS can hence be thought of as an assembly of pseudo ID quantum wires. 
The progress in this area has been the subject of many recent reviews*-^"^-* . In this paper 
we confine ourselves to the understanding of the optical absorption process in PS. 

There have been many attempts to obtain the band gap of PS from optical 
absorption measurements (9~i2)_ jj^ these measurements, the authors assume that the 
absorption coefficient (a) of PS satisfies the same relation as in 3D bulk c-silicon viz. 

Vahuj = A{nuj - Eg) (1) 



A plot of \J ah,uj vs h,uj is a straight line and the intercept on the energy axis 
gives the band gap Eg. The use of this equation for PS suffers from the following criticisms 

1) In general ,the absorption co-efficient depends on the joint density of states 
(JDOS) of the material. Eq.l is valid only when the density of states g(E) varies *^^^-' with 
energy near the band edges as g(E) = B \fE. If we treat PS as a ID material, it would be 
more appropriate to write g(E) =-^, where E is measured from the band edge. Clearly 
the JDOS in PS is expected to be different from that of bulk c-silicon. For a ID system 



having a single direct band gap Eg, it is easy to show that the JDOS is proportional to 
^ . In such a case, a plot of a{E) vs E should peak at Eg. On the otherhand, for an 



indirect band gap material ( assuming that momentum k is not conserved in this pseudo 
ID system during optical transition ), it is easy to show that the JDOS is a constant and 
independent of incident energy ( following reference (13) ) . Therefore , in either case , 



yafwo is not linearly related with huo at all. Clearly the use of eq. 1 to extract the band 
gap is not justified for PS . 

2) The nano structures that make up PS have a distrihutiorS^'^'^' of diameters 
(d). If the band gap is related to quantum size effects, then it is clear that the energy 
upshift Ai5 ( which can be written as AE = C/d"^ ) is different for each of the 
nanostructures that make up PS.( The value of X in the simplest approximation is 2 , 
further details of this are discussed in section III ). In PS, we are actually dealing with a 
heterogeneous system with a distribution of band gaps. Hence it may not be meaningful to 
visualize PS as having a single band gap and that this can be extracted from equation 1. 

In an attempt to address these questions we have investigated the optical 
absorption in PS .In this paper we clarify the absorption process in PS assuming that PS 
has a distribution of band gaps . Section II contains the experimental results. Here we 
report a non — destructive way of measuring the porosity of PS thin films using 
transmission measurements. Section HI outlines the model used for calculating optical 
absorption and section IV the results of our simulations . Finally we summarise our 
conclusions in section V . 

II. EXPERIMENTAL RESULTS 

PS is made by electrochemical etching of P (P'^) c-si in a 1:1 HP : Ethanol 
solution at a current density between 10-50 mA/cm^. Free standing samples are lifted off 
the substrate by increasing the current density to > 200 mA/cm^. The samples are then 
washed in water to get rid of HF , rinsed in n-pentane and finally dried in air on a glass 



substrate for transmission measurements . 

Transmission measurements are carried out using a CVI 240 Digikrom 
monochromator with a tungsten lamp as a source, silicon detector and SR 530 lock -in 
amplifier .For some samples we also use a Gary 1756 spectrometer. We have checked that 
internal multiple scattering does not dominate the absorption in our samples. This is done 
by making two PS samples of two different thickness {5fi and 10/i respectively) and 
ensuring that the transmission curves are the same for both samples at low absorption ^^^' . 

The absorption coefficient a (A) as a function of wavelength A is obtained from 
the normalised transmittance T(A) using 

^ {l-Rrexpi-a{\){l-P)t) 
^ ' 1 - R^exp{-2a{X){l - P)t) ^ ' 

The reflectivity (R) is obtained in the low a region using the relation R = 
(1-T)/(1+T). P is the porosity of the layer and t is its thickness. 

The porosity of PS is usually determined by gravimetric measurements which 
also destroy the sample. We now show that optical transmission experiment allows us to 
measure porosity in a nondestructive fashion . 

The normalized transmittance ( in eq.2) for a = reduces to 

(i-fi) i-n ,,, 

If n is the real part of the refractive index , then we can also write for the 
reflectivity ( for a = ) as 



{n-Vy 

(n + l) 



R-'^^^. (4) 



It follows from eq.4 that the refractive index {nps) of PS is 



l + Vfi „x 



We model PS as consisting of two media - air ( having a refractive index ( Uair 
= 1 )) and crystalline silicon ( ngi = 3.44 ) . If P is the porosity of a sample then we can 
write 

UairP + nsi{l - P) = nps (6) 

p ^ {nsi - nps) , . 

{nsi - 1) ^ ^ 

The porosity can now be obtained using equations 4 to 7. 

To check the validity of eq.7, we calculate the porosity ( using eq.7 ) and 
compare it with the corresponding values of P obtained gravimetrically. Figure 1 is a plot 
of the porosity determined from the transmission curves and the corresponding 
experimentally determined porosity by gravimetric measurements as reported in the 
literature'-^'^^-' . We see from figure 1 that the points nearly fit a straight line with slope 
one . We hence conclude that optical transmission measurements can reliably be used to 
obtain the porosity of porous silicon in a non — destructive fashion . 

Figure 2 shows the normalized transmission data on free standing PS films 
made from p type ( resistivity 4 ohm-cm ) si substrate and also from p'*' si( 0.01 ohm-cm ) 
substrate. We see from the figure that for the PS grown on p type si, the onset of low 
absorption occurs at a smaller wavelength than that for PS grown on p'*' si. This signifies 
that the p PS has a higher band gap than p+ PS and is in agreement with the reported 
results ^^'^\ 

Figure 3 is a plot of log(Q;) vs E at 300 K and at 100 K. The absorption at low 
temperature is reduced with respect to that at room temperature. We see that the two 
curves exhibit a rigid shift in agreement with other results '^^°-' . An explanation of this 
phenomena on the basis of our simulated results is given in section IV . 

We now attempt to understand the optical absorption process in PS with the 
help of some model calculations and simulations , assuming that PS has a distribution of 



band gaps . 

III. A MODEL FOR ABSORPTION IN PS HAVING A DISTRIBUTION OF 

BAND GAPS 

We model the nanostructures that make up PS to be ID parallelopipeds of 
square cross sections of side d and constant length L. Since d is typically 20 - 40 A, the 
band gap of the material is enhanced due to quantum size effects . We assume that the 
effective band gap ( lowest energy gap ) for a particular parallelopiped is the minimum 
separation between the conduction band and the valence band states for quantum numbers 
n^: = 1 and n^, = 1. We also assume that the absorption takes place in the bulk. 

The absorption coefficient a{E) for a particular value of incident energy E can 
be written as 

a{E) = f p\E - EG)F{EG)P{EG)dEG (8) 

J Ec,<E 

where the integration is done over all the band gaps (Eg) of the system below the incident 
energy E. In equation 8, p,j is the joint density of states, F(Eg') the oscillator strength for 
the optical transition and P(Eg) the distribution of band gaps Eg- in the material. We now 
discuss eq.8 in some detail. 

We have already pointed out in section I, if g(E) =B/yE then the JDOS is 
p^ = A{E — Eg)~^ where A is a constant. The value of W depends on whether PS is a 
direct or an indirect band gap material. For direct band gap material W = 0.5 and for 
indirect band gap material W = . 

F(EGr) is the oscillator strength governing the absorption. For indirect band gap 
material, the oscillator strength can be enhanced for small d'-^^-' . This happens because 
the overlap between the electron and hole wavefunctions in k space increases as a result of 
quantum confinement and contributes to an increase in the oscillator strength at small d 
(large Eg)^^^'^"^^ and can be written as F = i/d'^ where 7 = 5 to 6. This manifests itself as 
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the no phonon line which has been seen in the luminescence spectra of PS '^^^\ 

The third term in eq.8 viz. P{Eg), is the distribution of band gaps in PS. The 

band gap distribution is obtained by assuming a distribution of sizes for d and a relation 

governing the upshift in energy A E with the size d ( due to quantum confinement ). 
In our model we have considered two possible distribution of sizes. One is 

Gaussian and the other is Lognormal*^^^) distribution given by P'^(d) and P^(d) 

respectively. 



pG^d) = —^expi- ^"^ y' ) (9) 

where do is the mean size and the a the standard deviation. 



where mg = In(do), and a^ = ln(o"), do being the mean size and a the standard deviation. 

Electron microscopy measurements*-^'^'^'^-' suggest that for p PS, do is around 30 
A and a ^ 4 A. Figure 4 shows the distribution of sizes for Lognormal and Gaussian 
distributions for do = 30 A and a = 4 A. 

The energy up shift for the confinement AE can be written as 

AE = EG-Eg = C/d^ (11) 

where E^ is the crystalline Si fundamental indirect band gap ^1.17 eV and E^ is the 
increased band gap due to quantum confinement. 

The value of X is 2 in the usual effective mass approximation. However, for 
nanocrystals, the effective mass itself becomes size dependent and the exponent X has 
been reported to vary from 1.2 - 1.8 and becomes 2 at large size'-'^^'^^-' . For the purpose of 
our simulation we consider two cases 

a) X = 2, C = 486 eV/A^ (Ego = 1-71 eV for d = 30 A) (^0,22) _ 

and 
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b) X = 1.4, C = 126 eV/A^-^ (Ego = 2.25 eV for d = 30 A) (^i) . 

The interaction volume of the column structures with incident light is V(d) = L 
d^ . The resulting distribution R(d) = P(d) x V(d) is properly normalized to have equal 
number of dipole oscillators in each case. We now make a change of variable from d to Eg 
(Eg = Eg + C/d"^) to get the corresponding distribution P(Eg) of the band gaps. 

For a Gaussian distribution of sizes (eq.9) - the distribution of band gaps is 

P'^iEo) = Po r\Ea - EgY-^-^y^expi-Aii^^^y/^ - 1) ^i^G (12) 

Jo -&G — J^g 

where Pq is a normalization constant, A = 1/2 {do /a) and E;, is the maximum value of Eg- 
Following a similar procedure for the Lognormal distribution (eq.lO) we have 

PHEa) ^ R. f Vc - ^,).-v-.A..-e.H-B(||/^^ - D V^. (13) 

where Rq is the normalization constant , B = 0.5 (mo/o-i)^ . 

Fig. 5 shows the distribution of band gaps (P(Eg) when the energy upshift is 
given for X = 1.4 and X = 2. For the Lognormal case, we see that P(Eg) peaks close to 
the band gap of c-si. 

Equations 8, 12 and 13 are used to calculate a{E) for each incident energy E = 
hu. 

a{E) for a Gaussian distribution of sizes is given by 

a'^iE) = ao r\E - Ea)-^{Ea - i?,)(--^-^)/^exp[-A((f^^)V^) _ 1)^^^ (14) 
Jo HiG — ^9 

where a^ is a normalization constant . 

For a Lognormal distribution of sizes a{E) is 



a^E) = /^o/''\i5-i5c)-^(i^c-i^J«-^-^)/^V,^exp[-5(^^^-^;^-l)^]rfEG 

(15) 



where (3q is a normalization constant. 

Equations (14) and (15) are used to simulate various possible absorption 
processes . 



IV. RESULTS OF THE SIMULATION 

In this section, we present the results for the absorption co-efficient as a 
function of energy for different cases on the basis of the above model. The values of various 
exponents for different physical situations are summarized in Table I. 

(1) Direct band gap 

Figure 6 shows a{huj) vs huj for a direct band gap material where k is conserved 
in the optical transition and no size dependence of the oscillator strengths ( 7 = ) is 
considered (cases 1 and 3 and cases 5 and 7 in Table 1 ). To facilitate comparison between 
various cases , we have normalized the absorption to its peak value in each case. We see 
from the figure that the absorption tends to fall off at high energy ( in sharp contrast with 
the experimental data ). These results of the simulation can be easily understood as 
follows. 

For a ID system, the JDOS has a singularity at E = Eg and falls off for E > 
Eg- The dominant contribution (to each <y{E)) is hence only from a particular set of 
nanostructures with band gap Eg where E =Eg and all other Egs have almost no 
contribution. So the simulated a{E) vs E plot mimics P{Eg). We hence conclude from the 
simulation and the experimental data that PS is not a direct band gap material and we 
will not consider this case further. 

(2) Indirect band gap 

We now consider the case when momentum(k) is not conserved for the optical 
transition. Here we discuss two cases - size independent oscillator strengths (7 = 0) and 
size dependent oscillator strengths ( assuming 7 = 6). 



We first consider the situation where the oscillator strength has no size 
dependence. Figure 7 (curves (a) and (c)) is a plot of the calculated a vs h,uj for a Gaussian 
distribution of sizes (case 2 and 4 of table I ). Figure 8 (curves (a) and (c)) is a plot of the 
calculated a vs h,uj for a Lognormal distribution of sizes (case 6 and 8 of Table I). Also 
shown in the figures is a plot of the experimental curve for a p type PS. (The dip in a near 
2.2 eV in the experimental curve is an artifact related to a filter change). To facilitate 
comparison , the calculated a is scaled to match the experimental curve near 3 eV. We see 
from the plots 7(a) and 7(c), that a initially increases with energy and then saturates - at 
variance with the experimental results. We understand this result qualitatively as follows. 
In section I , we have already mentioned that the JDOS is a constant independent of 
energy for Tiuj < Eq and zero for huj > Eq- Hence for hu > Eq ,a given nanostructure 
contributes equally at all energies to a. However P(EGr) falls off as Eg increases. 
Nanostructures having a band gap {Eg) above a certain cut off in energy (in the high 
energy tail of the distribution) do not contribute significantly to a{huj) (as P(Eg') is very 
small). This gives rise to the saturation of a{E) in figures 7(a) and 7(c). This saturation is 
more obvious in figures 8(a) and 8(c), because P(Eg) falls off very rapidly for a Lognormal 
distribution of sizes. These results also do not agree with the experimental data. 

The situation is not remedied even if we include higher excited states (n^; = 2, 
Uy = 1; etc.) in the energy range considered in our calculation. Based on the above 
discussion, we see that for a to increase with hoo, it is necessary to increase the contribution 
to a from the nanostructures in the high energy region of a{E) .This is possible if the 
oscillator strength is enhanced for the nanostructures with large Eq ( smaller sizes ). 

We now consider the case when E{Eg) varies as f /d? . Figure 7 (curves (b) and 
(d) for case 2 and 4 in Table I )and figure 8 (curves (b) and (d) for case 6 and 8 Table I ) 
are plots of a{huj) vs hu for 7 = 6. For the Gaussian distribution that we have chosen, we 
see that the absorption still saturates although the onset of saturation moves to slightly 
higher energy. This is a measure of the contribution of the (high energy) tail of P{Eg) for 
the Gaussian case to the total absorption. For the Lognormal distibution, we see from 
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figure 8 (curves (b) and (d)) that a{huj) vs hu now qualitatively looks like the 
experimental curve. 

From a comparison of figures 7 and 8, we find that in order to get the 
absorption co-efficient to increase with incident energy , it is important that the 
contributions to a from the low energy regions of P{Eg) do not dominate the optical 
absorption at large hu. This is also satisfied if the distribution of band gaps P(Eg') is 
sufficiently broad. We illustrate this in Fig. 9 where we plot the calculated value for a for a 
Gaussian distribution of band gaps with the mean energy at 2.25 eV and a variance of 0.75 
eV. Hence we find that the nature of the band gap distribution P{Eg) plays a crucial role 
in determining the shape of the absorption spectrum . 

We conclude that to properly account for the optical absorption in PS we need 
to assume 

(1) PS is a pseudo ID indirect band gap semiconductor. 

(2) The oscillator strength is size dependent and is given by a power law f /d"^. 

(3) The distribution P{Eg) also plays an important role . Assuming a 
Lognormal distribution of the nanostructures diameters in porous silicon seems to 
naturally account for the absorption spectrum. 

All these considerations are in general agreement with results of luminescence 
experiments on PS .The low temperature luminescence of PS excited near the peak of the 
broad luminescence spectrum exhibits steps which have been associated with the 
characteristic TO mode of c-silicon ^'^^\ More recently'-^^-' well resolved peaks 
corresponding to the TO mode of Si have been reported for luminescence of PS suggesting 
that PS is an indirect band gap semiconductor. The observed dominant zero phonon 
component of luminescence constitutes direct evidence for the enhancement of the 
oscillator strength due to strong quantum confinement at smaller sizes (1^,23,24) ^ 

We now comment on the validity of the procedure for obtaining the band gap of 



PS from y/ahuj vs hu. Figure 10 is a simulated plot of \J ahuj vs h,uj ( case 8 in Table I ) 
for three different values of 7 and the experimental data on p type sample . First, we note 
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that the intercept of the hnear portion of each of the simulated curves is neither coincident 
with the peak position ( near 1.3 eV ) of the corresponding band gap distribution ( see 
figure 5 ) nor close to the mean band gap energy Eo = 2.25 ev. So we see that it is very 
difficult to relate this intercept with such physical quantities .Therefore , it is not obvious 
how to interpret the intercept with some kind 'effective band gap '. We also note that the 
intercept depends on the choice of 7. We hence believe that the intercept does not have 
the simple interpretation of a band gap as defined for a homogeneous system. We would 
like to point out that for such a heterogeneous system the gap can be interpreted only in 
an operational sense like E03 or E04 - depending on the value of the absorption coefficient 
(after correcting for the sample porosity). 

Finally we try to understand the temperature dependence of the band gap. To 
explain the rigid shift of the absorption curve with temperature (to lower a), we argue 
that Eg of each of the nanostructures follows a relation Eq (T) = Ego ~ iT, where Ego is 
the band gap at zero degree Kelvin . We assume that the change of the band gap with 
temperature is a consequence, like in bulk semiconductors, of the electron -phonon 
interaction and take 7 = 0.5 meV/K. In our model, changing temperature is equivalent to 
changing the value of C in the expression that determines the energy upshift AE. Figure 
11 is a plot of ln(a) vs huj for three different values of C. We see that the absorption curves 
are rigidly shifted vertically with respect to each other ^-^^^ as also seen in figure 2. We 
argue that this provides a natural explanation for the observed temperature dependence of 



a . 



V. CONCLUSIONS 

In this paper we have studied the absorption in free standing thin films of PS . 

We have calculated the absorption coefiicient as a function of energy assuming 
that PS can be modelled as a pseudo ID system having a distribution of band gaps . 
From our calculations, we show that it is neccessary to look at PS as an indirect band gap 
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material. To satisfactorily account for the absorption, we need to invoke that the 
oscillator strengths has a size dependence. We also show that a Lognormal distribution 
of diameters of the nanostructures that make up PS appears to account for the measured 
absorption spectrum of p PS . 

We argue that one can not use eq.l to obtain the band gap of such low 
dimensional materials like porous silicon . We also explain the temperature dependence 
of the absorption in PS on the basis of our model. 

We have also shown that porosity can be inferred non — destructively from a 
transmission measurement in the region of low absorption . 
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Gaussian distribution in d 


Lognormal distribution in d 


Case 


Value of X 


Nature of Transition 


Case 


Value of X 


Nature of Transition 


1. 


2.0 


Direct 


5. 


2.0 


Direct 


2. 


2.0 


Indirect 


6. 


2.0 


Indirect 


3. 


1.4 


Direct 


7. 


1.4 


Direct 


4. 


1.4 


Indirect 


8. 


1.4 


Indirect 



TABLE I. These physical situations we have considered in the model for the optical 
absorption in porous silicon .For direct transitions W = 0.5 and for indirect transitions W 
= . The value of 7 is taken as 6.0 for size dependent oscillator strengths and zero for the 
size independent case. 
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FIGURES 
FIG. 1. Porosity calculated ( using eq. 7 ) from the transmission curves'^'^^-* are plotted 
against the corresponding gravimetrically determined values of porosity . The points fall nearly on 
a straight line of slope one. (a) data from ref. 9 and (b) data from ref. 14 . 

FIG. 2. Normalized transmission spectrum for (a) p"*" and (b) p type porous silicon. The onset 
of low absorption is at smaller wavelength for p si than that of p^ si. 

FIG. 3. Temperature dependence of the absorption spectrum of p type porous silicon . (a) 
300 K and (b) 100 K . 

FIG. 4. Size distribution of the column widths . (a) For a Gaussian and (b) for a Lognormal 
distribution of widths, having do=30 Aand a= 4 A. 

FIG. 5. Band gap distribution for corresponding size distribution given in figure 4. (a) and 
(b) are for Gaussian distribution for X = 2 and 1.4 respectively, (c) and (d) are for the Lognormal 
distribution for X = 2 and 1.4 respectively. 

FIG. 6. Simulated a{E) vs E plot for direct transitions (7 = 0). Plots (a) and (b) are for 
the cases (1) and (3) in Table I . Plots (c) and (d) are for the cases (5) and (7) in Table I. 

FIG. 7. Simulated a{E) vs E plot for indirect transitions .(a) and (b) correspond to case (2) 
in table I for 7 = and 6 respectively .(c) and (d) correspond to case (4) in table I for 7=0 and 6 
respectively .(e) Measured a{E) vs E plot for p type PS. 

FIG. 8. Simulated a{E) vs E plot for indirect transitions .(a) and (b) correspond to case (6) 
in table I for 7=0 and 6 respectively .(c) and (d) correspond to case (8) in table I for 7=0 and 6 
respectively . (e) Measured a{E) vs E plot for p type PS. 

FIG. 9. (a)Simulated a{E) vs E plot for broad gaussian distribution of band gaps ( Eq = 2.25 
eV and cte = 0.75 eV ) with indirect transitions , X=1.4 and 7 = 6. (b) Measured a{E) vs E 
plot for p type PS. 
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FIG. 10. Simulated plots of {aTaoy/'^ vs fiuj for (a) 7 = 5 , (b) 7 = 5.5 , (c) 7 = 6.0 and (d) 
corresponding plot for the experimental data on p type PS . 

FIG. 11. Simulated \og{a{E)) vs E plots for various temperatures ( for case 8 with 7 = 6). 
(a) C = 126 = Co , (b) C = 1.25Co , (c) C = 1.5Co . Curves are rigidly shifted vertically from one 
another. 
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